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We investigate the nonlinear response of GaAs-based photonic crystal cavities at time scales 
which are much faster than the typical thermal relaxation rate in photonic devices. We demonstrate 
a strong interplay between thermal and carrier induced nonlinear effects. We have introduced a 
dynamical model entailing two thermal relaxation constants which is in very good agreement with 
experiments. These results will be very important for Photonic Crystal-based nonlinear devices 
intended to deal with practical high repetition rate optical signals. 
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I. INTRODUCTION 

With both a small modal volume and large quality 
factor, optical microcavities exhibit a greatly enhanced 
light-matter interaction and a strong non-linear optical 
response 0, [H, 0, H]- An emerging class of optical micro- 
cavities is based on air-clad two-dimensional (2D) Pho- 
tonic crystals (PCs). A Q-factor greater than 10^ was 
achieved with this technology [il, 12, which also al- 
lows small modal volumes (~ (A/n)^ or 10~^^m^. More- 
over, 2D PC technology is a planar technology which 
is particularly suited for the fabrication of photonic cir- 
cuits. Therefore, the move towards all-optical process- 
ing, long considered impractical, has changed dramati- 
cally and new possibilities have been opened [sj]. Im- 
pressive experimental demonstrations of low- energy ( fj ) 
optical bistability and all-optical switching [To, Till , 
wavelength conversion [12], optomechanical effects [l3] 
and dynamical control of the cavity lifetime [14] are re- 
cent noteworthy achievements. Most of these results 
come from a specific technology, i.e. silicon-based air- 
clad PCs where the nonlinear process involved is two- 
photon absorption (TPA) followed by plasma-induced 
and thermally-induced refractive index change. The op- 
tical power required is quite small (~ 100/ J) and can 
be very fast (70 ps by ion implant, as the carrier re- 
combination time result much shorter in these nanostruc- 
tured devices than in bulk silicon [l5|, [IB] • It has been 
predicted that, owing to strong light-matter interaction, 
nonlinear properties can be engineered by introducing 
nanoparticles [ij • 

The desire for an even faster response time motivates 
the research on alternative materials with a strong optical 
Kerr effect. Much progress has been made in processing 
chalcogenide crystals with high-Q PC microcavities have 
been demonstrated recently [18]. III-V semiconductors 
are also good candidates for optical switching as they 
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have two very attractive features, compared with silicon. 
First, the optimization of the Kerr effect with respect to 
TPA m, . Self phase modulation due to Kerr effect 
has been demonstrated in 2D PCs recently 0, An 
additional feature of III-V semiconductors is the possibil- 
ity of exploiting a strong nonlinear effect related to ab- 
sorption saturation in active structures, such as Quantum 
Wells (QWs) and Quantum Dots (QDs). Fast nonlinear 
dynamics [23|], leading to bistability [24] and excitabil- 
ity [25], has been observed in InP-based 2D PCs with 
InAsP QWs which are designed to operate with band- 
edge modes coupled to off-plane free space beams. Low 
power and very fast (2 ps) switching (15 ps for a com- 
plete on/off cycle) has been demonstrated with a Sym- 
metric Mach Zehnder-type all-optical switch made with 
InAs/AlGaAs quantum dots [26]. 

In this paper we focus on the dynamics of the pro- 
cesses initiated by TPA in GaAs PCs. Although, ideally, 
all optical switching requires the Kerr effect, in prac- 
tice TPA, followed by a carrier induced index change, 
is still a very attractive approach, because of the rela- 
tive simplicity of the technology. In particular, there is 
no active material, e.g. QDs, or phase matching condi- 
tion required. Moreover, compared to silicon, the TPA 
coefficient in GaAs is tenfold higher [23] and the non- 
radiative carrier lifetime in patterned structures can be 
very short (8 ps) [28]. A fast and strong nonlinear re- 
sponse is therefore expected in GaAs PC nanocavities. 
Very recently, we have demonstrated ultra-fast (6 ps re- 
covery time) and low power (^ 150 fJ) modulation in 
GaAs Photonic crystal cavities [29]. We have investigated 
optical bistability in high-Q {Q = 0.7 x 10^) PC cavi- 
ties on GaAs and reported an ultra-low threshold power 
(/iiy range) [sO, [IJ [32[. After the submission of this 
manuscript, memory operation have been demonstrated 
in a InGaAsP PC cavity[33] 

In these experiments, only the slow {/as) regime has 
been explored, which is dominated by the thermally in- 
duced index change. This is also due to the high ther- 
mal resistance of the membrane structures compared to, 
for instance, micro-disks [34] or membranes bonded to a 
Si02 cladding [2j] and the lower thermal conductivity of 
GaAs with respect to silicon. 
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In this paper we investigate the response of PC micro- 
cavities at a much faster modulation rate (up to 50 MHz), 
which is weh beyond the typical thermal relaxation time 
of photonic devices. Under these conditions, the analysis 
of nonlinear responses, such as bistability, cannot be ex- 
plained with static models. First of all, moving to faster 
time-scales modifies the relative infiuence of thermal and 
carrier plasma effects. While at very fast time scales (ps) 
the dynamics tend to be controlled by the carrier lifetime, 
we will show that in the range between 1 - 100 ns the dy- 
namics results from the interplay of fast thermal effects 
and carrier plasma index shift. We introduce a model 
that incorporates two thermal relaxation constants. This 
is necessary to explain this dynamics and the fact that, 
despite the high thermal resistance of PCs microcavities, 
there are thermal effects which develop in less than 10 
ns. Section II is devoted to experiments made on a PC 
microcavity coupled to a waveguide. In section III we 
will introduce our multi-scale model. Discussion of the 
results forms the body of section IV. 



II. EXPERIMENTAL SETUP AND SAMPLE 
DESCRIPTION 

The photonic crystal structure studied here is the well- 
known optimized three missing holes (L3) PC microcav- 
ity [35[ in an air slab structure (thickness is 265 nm) 
based on a triangular lattice (period a = 400nm) of holes 
with radius r = 0.24a. The holes at the cavity edge were 
shifted by s = 0.15 a. The cavity is side-coupled (the 
spacing is 3 rows) to a 1 mm long line-defect waveguide 
along the TK direction, with width W = 1.05 a/3 a (see 
figure [1]). The fabrication process and the detailed linear 
characterization of similar structures are described else- 
where 0. The loaded and intrinsic Q-factors, 7,000 and 
30,000, respectively, are estimated from measurements 
using the procedure discussed in a previous paper [37l]. 
The cavity resonant wavelength is 1567 nm. The charac- 
terisation setup consist of a tunable external cavity semi- 
conductor laser (Tunics) with a relative accuracy of 1pm 
and narrow linewidth (<< 1 MHz), which is amplified by 
an Erbium Doped Fiber Amplifier (EDFA). The polarisa- 
tion is set to TE, or electric field in the plane of the slab. 
Coupling to the PC waveguide is obtained through mi- 
croscope objectives (Zeiss, N.A. 0.95) and microposition- 
ers. The transmitted signal is detected with an InGaAs 
photodiode (bandwidth = 2 GHz), amplified (40dB, tran- 
simpedance, 2GHz bandwidth) and monitored by an os- 
cilloscope (Tektronix 12.5 GS/s). 

Two distinct regimes of modulation are considered: 1) 
sinusoidal modulation (frequency from few kHz to 50 
MHz) and 2) low duty cycle (20 ns) square pulses. The 
cavity is observed from the top with a microscope (long 
working distance objective Zeiss, N.A. 0.4) and an In- 
frared InGaAs camera (Xenics). When the cavity is on 
resonance a bright spot is observed. 
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FIG. 1: (Color online) Top: SEM image of the sample, also 
showing a line-defect photonic crystal waveguide with side 
coupled L3 cavity. Bottom: linear transmission near reso- 
nance showing that the cavity is over-coupled. The linewidth 
of the cavity is larger than the transmission dip because of 
interplay with the Fabry Perot fringes. The left scale is the 
waveguide transmission and the right scale is the signal de- 
tected by the IR camera. 



A. Sinusoidal modulation 

At a low modulation rate (kHz) it has been shown by 
several groups [ol, ITol, ITll, [31I, that bistability occurs at 
a fairly low optical power (/iVF range, coupled into the 
waveguide). The dominant nonlinear effect is thermal 
induced index change because of heating resulting from 
carriers generated by TPA. Since this effect leads to a red- 
shift of the cavity frequency, bistability is observed when 
the initial detuning AA|t=o = Xiaser - Kavity\t=o) is pos- 
itive, where Xcavity\t=o is the cold cavity resonant wave- 
length. This is shown in fig. [2l Typically, the detuning 
is set between V3/2 (the theoretical minimum, although 
it has been shown that the Fabry- Perot resonances in 
the waveguide affect this value [ll|) and several times 
the cavity linewidth, estimated to be AXfwhm = 220 
pm in our case. If the modulation frequency is increased 
while the detuning and the average signal power are kept 
constant, the bistable behavior disappears. For instance, 
with a detuning ^ 330pm ^ l.bAXpwHM and a mod- 
ulation frequency of 15 kHz, bistability is observed in this 
sample at SO jaW. The power coupled into the waveguide 
is estimated considering the amplified laser power and 
the input objective coupling factor (here ~ -10 dB) as in 
our previous work (3l|. 

When the modulation frequency is increased to 1 MHz, 
the bistable effect and all traces of nonlinear distortion 
of the transmitted signal disappear. This could be ex- 
plained in terms of the slow response of the thermo- 
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lifetime of carriers in GaAs compared to silicon dras- 
tically reduces the strength of this effect, compared to 
thermo-optic effects. 
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FIG. 2: (Color online) Low modulation frequency experiment. 
Oscilloscope traces (a) showing the typical bistable response 
of a PC microcavity at low modulation frequency (15 kHz). 
The average input power in the PC waveguide is about 80 /xW. 
The bistable cycle is also sketched and the transitions 1-4 are 
also marked. Detunings are +0.4, +1.5 and +1.55AAFWi^M 
respectively. Scheme representing a cavity with a negative 
detuning b) and going to resonance due to nonlinearly induced 
red-shift c). 



optical effect. An important question arises whether the 
carrier plasma effect, much faster but weaker, will then 
take over. To answer that, the signal power and the 
modulation frequency were increased further and the de- 
tuning (still positive) was swept continuously from zero 
to about +700 pm, with the laser always on. This de- 
tuning value is larger than in low frequency experiments. 
A strong nonlinear distortion of the transmitted signal 
is then observed. The onset of the nonlinear behav- 
ior is simultaneous to the observation of a bright spot 
through the IR camera, thus indicating that the reso- 
nant frequency of the cavity is red-shifted under a sinu- 
soidal modulation pump. In particular, fig. [3] shows a 
typical output signal when the modulation frequency is 
increased to 50 MHz and the average power coupled in 
the waveguide is about 2.5 mW. The thermo-optic effect, 
driving the bistability at low modulation frequency, in- 
tegrates over time. This explains ultra- low power levels. 
At higher modulation frequencies the power required to 
observe nonlinear effects increases. The nonlinear behav- 
ior is observed with a positive detuning from 600 pm and 
670 pm (^ 3AAfw^//m)- We don't believe that the 
observed nonlinear distortion corresponds to a bistable 
behavior. Indeed, a fast nonlinear mechanism (plasma 
induced index change) should be dominate all other non- 
linearities. That is probably not the case, as the short 
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FIG. 3: (Color online) High modulation frequency experi- 
ment: oscilloscope traces of the transmitted signal with si- 
nusoidal modulation (50 MHz) for different detunings. The 
average power is about 2.5 mW. 



The nonlinear distortion observed in figure [3] is quite 
general and can be reproduced at different modulation 
frequencies in the 10 - 100 MHz range. The nonlinear 
transitions are not sharp, compared to the modulation 
period (response time is a few ns). This time is however 
much faster than what could be attributed to a ther- 
mal effect. Additionally, the period, T = 20 ns, is much 
shorter than what is estimated to be the thermal recovery 
time in air slab PC cavities [9]. We claim that the faster 
carrier plasma effect plays a key role at the ns time scale, 
while the thermal effect accumulates to strongly redshift 
the cavity. 

This mechanism is represented in Fig. [H Initially (a), 
the cavity is at room temperature and laser is detuned 
negatively (AA < 0). The laser wavelength is increased 
continuously (b). The cavity is blueshifted (plasma ef- 
fect) very fast, but heating quickly dominates and keeps 
the cavity redshifted until the end of the modulation cy- 
cle. Therefore, the cavity will be pushed towards longer 
wavelengths as long as the heat generated by TPA over 
each cycle compensates the thermal flow across the mem- 
brane. Over each modulation cycle, the cavity resonant 
wavelength will oscillate about an equilibrium point set 
by this thermal balance alternating blue-shift and fast 
heating (c). The strongest distortion of the sinusoidal 
modulation is observed nearest to the largest detuning 
which can be sustained. In our case, we achieved detun- 
ings of up to 700 pm, which is related to a local tem- 
perature increase of 6 K through the thermo-optical co- 
efficient (table. H]). Better physical insight into the two 
competing effects (thermal and free carrier), is gained by 
reducing the average input power, e.g considering much 
smaller duty cycles. 
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FIG. 4: (Color online) Scheme representing the interplay of 
thermal and carrier plasma effect and how they are related 
to the detuning. The system is at rest (cavity at its cold 
frequency) and the laser is switched on (a). The laser wave- 
length is increased but heating keeps the cavity red detuned 
(b). When the cavity temperature is high enough, the mean 
resonant wavelength get closer to the laser and the nonlinear 
distortion is apparent (c). 



B. Low duty cycle excitation 

The cavity is excited with low duty cycle, square pulses 
(duration At = 20 ns, period T=2 /is), so that the cav- 
ity has time to cool down before the arrival of the next 
pulse. We refer to this case as a "non-heating pulse" 
because the average temperature of the cavity remains 
close to room temperature. Fig. [5] reports the input and 
the transmitted pulse depending at various wavelength 
detunings. The input power in the pulse is estimated to 
be on the order of 10 mW. A more accurate determina- 
tion of this value was complicated by the very low filling 
factor and the use of the EDFA, which is not designed 
to operate in these conditions. This lead to some fluc- 
tuation in the level of the amplified pulse. Furthermore, 
very low duty cycle means that the the power contained 
in the pulse is a fraction of the amplified spontaneous 
emission(ASE), which means that it is difficult to deter- 
mine the exact value of the measurement of the average 
power. Let us first consider a negative detuning: e.g. 
A A ^ —0.5AXfwhm (fig- W^)- The leading edge of the 
output pulse rises with almost linear slope and is about 
5 ns long. For zero or moderate positive detuning (e.g. 
A A between and ^ A^fwhm)^ the output pulse typ- 
ically has a two step shape with the leading part of the 
pulse having a lower level (fig. [Sb and c). The transi- 
tion between the two levels is sharp (a few ns) and moves 
from the leading edge to the trailing edge of the pulse as 
the detuning is increased. If detuning is increased fur- 
ther, the pulse takes on an almost triangular shape with 



a sharp leading edge and a long trailing edge (figure [5] d) 
before losing any signature of nonlinear distortion. 
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FIG. 5: (Color online) Experiment with "non heating" pulses. 
Input (dashed) and output (solid) pulses as a function of the 
detuning, (a) AA — —0.5AXfwhm , (b) -\-0.7AXfwhm , (c) 
-\-0.9AXfwhm and (d) -\-1.5AXfwhm . The pulse period is 
2/js. 

Interestingly, the initial leading edge of the output 
pulse is very similar in the four cases considered in fig. 
O i.e. we don't see low transmittance associated to the 
on-resonance case, as long as the power is high. We see 
this behavior at much lower power only. We believe that 
this is due to a fast carrier effect which detunes the cavity 
faster than our detection apparatus is able to measure (1 
ns). 

It is also important to note that no signature of a 
plasma-induced bistable state is observed when the de- 
tuning is negative. We think that this is because of fast 
heating, which dominates the plasma effect before the 
pulse is extinguished. These features are well explained 
by our model. 



III. NONLINEAR DYNAMICAL MODEL 

The coupled-mode model developed here is based on 
previous optical microcavity literature. Particular struc- 
tures that have been investigated include: microdisks js^ 
and photonic crystal cavities [H, [10, |4l| . The dynamical 
variables are: the optical energy inside the cavity |a(t)p, 
the free carrier density N{t) and the cavity and mem- 
brane temperature T{t). The model considers a single 
mode [47] cavity (resonance is at ujq and unloaded qual- 
ity factor Qo), which is side-coupled to a waveguide (the 
loaded Q-factor is Q). The modulation of the transmis- 
sion observed in figH] is related to the Fabry-Perot reso- 
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nance due to the finite reflectivity at the waveguide end 
facets. This spatially extended resonance has no impact 
on the nonlinear response (other than merely modulating 
spectrally the coupling into the cavity), since the associ- 
ated fleld intensity is orders of magniture smaller than in 
the cavity. The fleld in the cavity follows the equation: 



da{t) 
dt 



{iuj — iujL 



W) + \^Pin{t) (1) 



where a;^ is the laser frequency, Tc = 00{){Q~^ — Qq^) 
gives the cavity to waveguide coupling strength, Ttot is 
the inverse (instantaneous) cavity lifetime, Pinif) is the 
power in the waveguide and uo = uoo + Acj y is the instan- 
taneous cavity frequency. Following |3l, li^ , the nonlin- 
ear change of the cavity frequency is given by: 



NL 



-An : 



712/ C 



rieff uVk e 



\a{tt 



-AT(t) + -Nm 



(2) 



Here n is the refractive index of the bulk material and 
rieff the effective refractive index in the cavity, i.e.: 



eff 



J n{r)^ \E{r)\^ dV/ J \E{r)\UV, and n2i is the 



Kerr coefficient and Vke 
defined as: 



the Kerr nonlinear volume. 



Kerr 



Jn2i{f)/n2i\E{f)\^n{r)^dV 
{Jn{rr\E{r)\'dVr 



(3) 



The refractive index change due to the plasma effect is 
dn/dN = — cj^/2ncj^A^, with cOp = e^A^/eom* the plasma 
frequency. In principle, holes and electrons both con- 
tribute to this effect, however, given the much smaller 
effective mass of electrons in GaAs, the contribution of 
holes is negligible. The inverse instantaneous photon life- 
time is: 



tot 



TP A 



^FCA 



(4) 



the first term is the inverse linear cavity lifetime, Ttpa 
and Tfca are the contributions from two photon (TPA) 
and free carrier absorption (FCA), respectively. The 
TPA term is: Ttpa = p2C^ /n'^Ht)\^/VTPA with p2 
representing the TPA coefficient in units of m/W and 
Vtpa = Vxerr the nonlinear effective volume [40], while 
the free carrier absorption is proportional to the com- 
bined free carrier density TpcA = {cre-\-(Jh)N{t)c/n. The 
evolution of carrier density follows the rate equation: 



dN{t) _ N{t) c^/n^P2 



1 



dt 



TN 



Vtpa 2huj Vcc 



(5) 



Here h is the reduced Planck's constant, tn is the effec- 
tive carrier lifetime, and Vcar is the volume in which the 
carriers spread and recombineji^]). In this approxima- 
tion, the population decay is dominated by recombina- 
tion at the surface, due to the large surface to volume 



ratio typical of photonic crystal strucures. Carriers are 
assumed to spread and distribute homogeneously within 
the carrier volume Kar, which is assumed to correspond 
to the region of the membrane delimited by the holes 
around the cavity. This approximation is rough, but it is 
a reasonable choice. 
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FIG. 6: (Color online) Modelling heat diffusion in the PC 
membrane. Local cavity temperature Tc vs. time calculated 
with a 2D thermal diffusion equation (dots) and fitting (solid 
thick and thin lines) with the model eq. [6] Inset: enlarged 
view of the short and the long term behaviour. Inset right: fit- 
ted rth,m = 200ns with a single exponential (dashed line). Fit 
with 4 unconstrained parameters (thin line). Constrained fit 
with Tth,m = 200ns rth,c = 8.4 ns (thick solid line, main 
plot and inset left). Thermal capacitances for the cavity and 
the membrane are 0.43 10"^^M^/i^^ and 3.4 10"^^l^/i^^ respec- 
tively. 

Because of the high thermal resistance in membrane 
PCs, thermal effects are very important and must be 
modelled appropriately. This is, for instance, the scope 
of ref.[42[. In contrast with existing literature, we intro- 
duce a more complicated model for the thermal response. 
We assume that the heat generated in the cavity has 
not yet diffused to the border of the membrane at the 
time scale of interest (1 ns). Therefore, a very impor- 
tant physical quantity is the thermal capacitance of the 
cavity Cth,c = CypVth,c^ associated with a small region 
of the membrane, with volume Vth,c roughly correspond- 
ing to the cavity volume. Heating and the spreading of 
the heat are modelled by solving the two dimensional 
heat diffusion equation. The heat is generated for some 
time (10 ns) inside the cavity and then the system is al- 
lowed to cool down. We assume that the radiative and 
convective contributions are negligible; therefore, since 
the membrane is suspended in air, all the heat has to 
ffow through it. The result of this calculation is shown 
in fig. [6l The cavity temperature Tc increases with a 
rate that is governed by Cth,c and Tc decreases at two 
exponential time scales. The fast time scale Tth,c is asso- 
ciated to a relatively fast transfer of heat from the cavity, 
where it is generated, to the neighbouring region in the 
membrane. We associate a second thermal capacitance 
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TABLE I: Physical paramenters used in 


the dynamical model. 


Parameter 


Symbol 


Value 


Ref. 


TPA coefficient 


f32{cm/GW) 


10.2 


m 


Kerr coefficient 


n2i{c7n'/W) 


1.610"^^ 


\27] 


Loaded Q 


Q 


7000 


meas. 


Intrinsic Q 


Qo 


30000 


est. 


Modal Volume 


^mod 


0.66(A/no)^ 


calc. 


TPA Volume 


Vtpa 


3.13(A/no)^ 


calc. 


Carrier Volume 


Vcar 


7 Vmod 




Thermo-optic coeff. 


dn/dT (K-^) 


2.48 10"^ 


[43] 


Therm, eff. vol.(c) 


yth,c 


3.8 * Vmod 


calc. 


Therm, eff. vol. (m) 


Vth,m 


Q1 ^ T/ 

* Vmod 


calc. 


Specific Heat 


Cvp {W/Km- 


1.84 10^ 


[43J 


Carrier lifetime 


TN 


8 ps 


[28, 29] 


Therm, relax, time (c) 


Tth,c 


8.5 ns 


calc. 


Therm, relax, time (m) 


'Tth,m 


200 ns 


calc. 


Therm, resistance 




7.510^ K/W 


calc. 


FCA cross section 


ae,h (lO-^^m^ 


9.3(^)- 


[44] 



Cth,m to the whole membrane. The second time scale 
nh,m takes into account the spreading of the heat over 
the rest of the membrane (with temperature T^) to the 
bulk semiconductor structure. This picture is substan- 
tially different from previous models with a single time 
constant r = Cth,cRth that is obtained from the thermal 
resistance Rth, defined as = Rth^T (ratio between 
the increase of the temperature over the heating power 
at steady state). We will show that these two time scales 
play a crucial role in understanding the system response 
under sinusoidal and single pulse excitation. The corre- 
sponding model entails therefore two auxiliary equations: 



dT,{t) ^ T,{t)-Tm{t) , 2hLj V,, 
dt 

dTm{t) 



Tth,c 
Tm(t) 



-N{t) 



To 



Tat Cth,cVth 
Ti{t) — T2 Cth,c 



dt 



nh,i 



th,m 



(6) 
(7) 



The thermal effective volume of the cavity and of the 
membrane Vth.c Vth,m and their thermal relaxation times 
^th,c^ ^th,m are obtained by fitting the solution of the 
heat diffusion equation with the solution of the two time 
constant model (eq. [6]). If we fit all four parameters (two 
time constants and two capacitances) the fit turns out to 
be very good in the first 50 ns, yet tends to underestimate 
the long time constant and therefore to underestimate 
the total thermal resistance. To prevent this, we first fit 
a simple exponential to the long term behaviour, which 
gives a time constant of approximately 200ns. Then we fit 
the complete curve by constraining rth,m = 200ns. The 
result is shown in fig. [6] and confirms that this model 
is well adapted to describe the thermal behaviour of PC 
microcavities. The result of this calculation is given in 
table H 



IV. DISCUSSION 

Simulations have been carried out with parameters 
given in Table [H All of them are well known or can be cal- 
culated or measured with reasonable accuracy. The car- 
rier volume Vcar is the exception. The impact of carrier 
diffusion has been investigated theoretically in a recent 
paper [45], providing some hints to explain a fast recov- 
ery time in Silicon PCs. The carrier lifetime in patterned 
GaAs (e.g. 2D PCs) is however much shorter than in Sil- 
icon. Very recently we estimated it to be about 6 ps in 
our GaAs cavities [2y, which is consistent with previous 
(8 ps) estimates for GaAs PC structures [28]. 

In the limit where the dynamic is much slower than 
Tat, which the case considered in this work, tn and Vcar 
play the same role and what matters here is the ratio 
^N/Vcar- To show that, let us consider the density of the 
generated carriers : 



N{t) 



Vtpa 2fkd Vcar 



(8) 



The blue-shift of the cavity resonance due to the gen- 
erated carriers is therefore proportional to the instan- 
taneous power absorbed. Pa. The red shift induced by 
heating is proportional to the absorbed energy and to 
the inverse of the thermal capacitance. Then, the rela- 
tive strength of thermal and carrier effects depends on 
the ratio tn /Vcar and on the thermal capacitance. The 
dynamical behaviour of the system investigated here is 
basically determined by these two quantities. 

Modelling was performed by varying only one param- 
eter, the carrier volume, and adjusting the input power 
and the detuning around reasonable values. All other 
values are well known or calculated precisely. 

Fig. [Tlreports the simulated responses of a square pulse 
at various detunings. When the initial detuning is nega- 
tive (e.g. AA(0) = —120pm, fig. [3 a and b), the carrier 
effect is responsible for the drop of the leading edge of 
the transmitted pulse. The signal coupled in the cavity 
is strong enough to cause considerable generation of car- 
riers and the resulting blue-shift tends to tune the cavity 
into resonance so that transmission is reduced (the cav- 
ity is side-coupled). Within a ns time scale, the heating 
following carrier recombination red-shifts and therefore 
detunes the cavity again. Such fast heating arises from 
the small thermal capacitance of the cavity. After the 
pulse is extinguished, the cavity is red-shifted by about 
100 pm with respect to the initial position and recovers 
its initial state in a fraction of jas. Let us now consider 
the case of zero detuning (fig. [7t and d). The carrier 
induced index change produces an instantaneous (with 
respect to the time scale considered here) blue shift that 
adds to the initial detuning. Heating, following absorp- 
tion, tunes the cavity back and transmission decreases 
once again. Thus more carriers are generated and the 
plasma effect tends to oppose heating but, in the end, 
the heating dominates. At these modulation rates ther- 
mal relaxation is not sufficiently fast to dissipate heat 
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time (ns) time (ns) 

FIG. 7: (Color online) Simulation of the response (solid) to a 
single pulse (dashed), depending on the detuning AA(0): -120 
pm (a,b), pm (c,d), 80 pm (e,f) and 150 pm (e,f). Right: 
corresponding instantaneous frequency shift AXnl of the cav- 
ity resonance (solid line) and carrier plasma (thin-dashed) and 
thermo-optical contributions (dotted). The laser frequency is 
also shown (thin solid). Pin = 1.2 mW, 0.9 mW, 1.2 mW and 
1.0 mW respectively. 



from the cavity after each duty cycle. Thus heat accu- 
mulates over time, causing a net red-shift of the cavity 
resonance. When the instantaneous detuning AXnl be- 
comes positive, carrier generation starts to decrease and 
so does the carrier induced shift. This initiates a positive 
feedback that quickly detunes (red) the cavity and makes 
the carrier density drop very fast. The delay between the 
pulse leading edge and this transition is clearly related to 
the amount of the initial detuning (fig. [7^ and f). When 
the initial detuning is positive and large enough (fig. [7^ 
and h), then heating tends to tune the cavity into reso- 
nance, thus reducing the transmission but also generat- 
ing carriers. The plasma effect will partially oppose the 
red-shift, thus explaining an almost linear change with 
time. The sharp step observed in fig. [5)3,c can only 
be explained by the thermal effect overtaking the carrier 
plasma index change. This step cannot be reproduced 
for any choice of free parameters if the plasma induced 
index change is suppressed in the model. This is shown 
in figure [HI In panel a and b we ran the same simulation 
as in fig. [7^, but we have suppressed the contribution of 
carriers to the frequency shift of the cavity. In this case 
the transmission increases with a linear, but finite, slope 
after the dip, independent of the pulse power. This is 
not what we see experimentally. Conversely, when the 
plasma effect is included, the sharp transition is repro- 




10 20 30 10 20 30 

time (ns) time (ns) 

FIG. 8: (Color online) Simulated response when the contri- 
bution to the plasma shift from carrier plasma is suppressed 
(a) and full model (c). The excitation power is 1.2 mW (thick 
line). The case with Pin =1.1 mW and 1.3 mW are also plot 
(thin lines). The detuning AA(0) is 80 pm. The correspond- 
ing instantaneous frequency shifts AXnl is also plot (b,d) for 
the case Pin =1.2 mW. 



duced correctly. There is a physical reason for that: the 
observed steep step requires a fast index change with, 
say, a negative sign, and a slower index change, with 
opposite sign. A simple explanation is that at first the 
plasma effect dominates (over a very short time, a few 
ns) before thermal heating takes over. If a plasma effect 
is ruled out, there is no way to explain the experimental 
results. A Kerr effect cannot be considered because it 
has the same sign as the thermal effect. After the sub- 
mission of this manuscript, critical slowing down (CSD), 
implying a transition from a low (off-resonance) to a high 
(on-resonance) transmission state, has been reported in 
InP-based PC cavities [46]. This effect also results into a 
step in the pulse response. We believe however that this 
is a different mechanism than what we describe here. The 
reason is that our system is side-coupled, thus CSD would 
manifest as a transition from a high (off-resonance) to a 
low (on resonance) transmission state. We did not ob- 
serve that. A possible explanation is the carrier lifetime 
in GaAs PCs being mush shorter than in Silicon or InP 
PCs. This implies that the plasma dispersion effect is 
much weaker than the thermal induced index change. 

We found experimentally that the response of the cav- 
ity is very sensitive to the power of the pulse, in particular 
in the case reported in [5^. This behavior is well repro- 
duced by theory. This means that input power must be 
set precisely in simulation and it cannot be the same for 
all the detunings considered here. Indeed, experimen- 
tally the peak pulse power is not constant since, a) the 
output power from the amplifier showed a marked depen- 
dence on slight change in wavelength, b) the waveguide 
transmission is modulated by a strong Fabry-Perot ef- 
fect (fig. [1]). The fact that the detuning in experiments 
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FIG. 9: (Color online) Simulation of the response of the cavity 
to a sinusoidal excitation for the following detunings (from 
top to bottom): 440 pm (a,d), 490 pm (b,e) and 500 pm 
(c,f). Input peak power is 1 mW. The transmitted pulses are 
shown on the top (a,b,c), the instantaneous frequency shift 
(solid think curve) and the thermo-optical contribution alone 
are shown on the bottom (d,e,f). 



and in the theory are not the same can be understood by 
accepting some error in the measurement of the detun- 
ing. The reason is that the laser was directly modulated 
which induced some deviation from the nominal setting 
of the laser on the order of 100 pm. 

This relatively fast (few ns rise-time) nonlinear effect 
reported in fig. [5] is well understood through modelling. 
The complexity of the model is justified by the number of 
different time scales present. An important point of this 
paper is to demonstrate that the same model is able to 
explain also the response to a sinusoidal excitation. The 
response is calculated for a time long enough to ensure 
that we reach the steady-state regime. The simulations 
shown in fig. [9] are carried out in order to mimic exper- 
imental conditions, that is, the excitation wavelength is 
detuned negatively with respect to the cold cavity fre- 
quency. As the excitation is turned on, the wavelength 
is increased adiabatically with respect to the modulation 
period. Indeed it is found that the cavity stays red de- 
tuned, and is pushed towards longer wavelengths, until 
the excitation power and the heat generated is enough to 
sustain the detuning. 

The simulations clearly support the experimental re- 
sults shown in figure [3] and interpretation in fig. [H In- 
deed, if we now compare the measured and calculated 
oscilloscope traces we conclude that the dynamics are ac- 
curately reproduced, as the detunings used in modelling 
correspond to what is measured experimentally. 

In particular, the response is very well understood 
by looking at the instantaneous frequency of the cav- 
ity, which crosses the laser frequency at different points 



depending on the detuning. If the detuning AA is too 
small, the cavity wavelength is always higher than the 
laser frequency, as the instantaneous cavity nonlinear 
shift AAatl = Kavity{t) - Xcavity{0) duc to heating is 
much larger over the whole modulation cycle. When the 
detuning A A increases, AAatl also increases on average, 
as the resonance is closer to the laser frequency, but less 
than AA does. Thus, there exists a combination of pulse 
power, modulation frequency and detuning such that A A 
and AAatl are close one to each other. If A A is increased 
further, heating is not enough to keep AAatl close to it 
and nonlinear distortion is lost. The role of the plasma 
induced index change and the interplay with the thermo- 
optic effect is also clear in figure [3^ b,d,f,). 

It is interesting to note that the plasma effect van- 
ishes when t < Ans and t > 16ns in figs. [9)3, d and f, 
although minima of the sinusoidal excitation are at t=0 
and t=20 ns (figs. [9^,c,e). This important point is un- 
derstood when considering that while the plasma effect 
is almost instantaneous (within the time scale relevant 
to this work), the thermo-optic effect follows the time 
integral of the optical power. At time t = the carrier 
population is negligible and the cavity is cooling down, 
and therefore is approaching the laser wavelength. As the 
power inside the cavity increases, carriers are generated 
with a rate proportional to the power squared. A carrier 
induced index change builds up very fast and produces 
the difference between the thin and the thick curves. The 
power absorbed produces heating. Instantaneous heat- 
ing is the time integral of instantaneous absorbed power 
and therefore builds up with a delay with respect to the 
plasma effect. The thermal effect then dominates the 
plasma effect. As soon as the resonance is crossed again 
and the detuning becomes comparable to half the cavity 
linewidth (100 pm), nonlinear absorption decreases very 
quickly and the carrier population declines accordingly. 
Therefore the carrier effect disappears. The cavity again 
starts to cool down and the cycle begins anew. 



V. CONCLUSION 

Switching based on TPA has a great potential in GaAs 
based structures, as it benefits from the strong nonlinear 
absorption coefficient (32 and sub-ns carrier lifetimes typi- 
cal of this material. However, most practical applications 
will require a high repetition rate (i.e. high duty cycle), 
therefore thermal effects are unavoidable in such small 
structures. In particular, heating is very strong in pho- 
tonic crystals, because of the poor thermal conductivity 
of PC membranes. We have experimentally investigated 
this regime and confirmed results from other groups evi- 
dencing that the thermal dynamics in PC membrane cav- 
ities is much faster than in macroscopic photonic devices 
due to the small modal volume of the structure. More 
particularly, we have explored a regime in which a sharp 
(2 ns) transition result from a positive feedback between 
thermal and carrier induced refractive index changes. 
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The peak power required to see these effects is at the 
mW level. We have also explained why it is not possible 
to observe bistability in GaAs cavities when the modu- 
lation rate is on the order of 10 - 100 MHz. Based 
on these experimental results, we have introduced a dy- 
namical model for PC microcavities which includes two 
thermal relaxation constants in order to account for heat 
diffusion from the cavity to the neighbouring membrane. 
We show that this is crucial in order to understand the 
dynamical response of the cavity and that this model re- 
produces the experimental results quite well. The ability 



to properly account for very different time scales is cru- 
cial for modelling patterning effects (i.e. long term or 
memory effects) which result not only from the carrier 
dynamics but also on fluctuations of the average power 
which induce changes in the local temperature of the cav- 
ity. This is crucial to design PC based all-optical switches 
for practical optical signals at high repetition rate (e.g. 
10 GHz). 

One of the authors (C. Husko) thanks the Fulbright 
Grant for financial support. 
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